Memory-effect in glasses at low temperatures 
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The dielectric constant of amorphous solids at low temperatures is governed by the dynamics of 
tunneling systems, small groups of atoms which tunnel between quasi equivalent potential minima. 
Recent experiments showed that at temperatures below 20 mK various glasses exhibit memory 
for a previously applied electric bias field. A first sweep of an electric bias field may prepare resonant 
pairs of tunneling systems, which are temporarily formed during the sweep, in metastable states. In 
subsequent sweeps the same resonant pairs thus significantly contribute to the dielectric constant, 
leading to a higher dielectric constant. We investigate the dynamics of resonant pairs during a bias 
field sweep yielding a qualitative explanation of the memory effect. 

PACS numbers: 61.43.Fs, 77.22.-d, 75.10.Jm, 05.70.LN 
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Glasses at low temperatures exhibit properties which 
differ considerably from those of their crystalline coun- 
terparts § . They are thought to be caused by tunneling 
systems (TSs) stemming from the tunneling motion of 
small atomic entities in double- well potentials. The en- 
ergy splitting E of these TSs is given by E 2 = A 2 + Aq, 
where A is the asymmetry energy due to the difference in 
the depth of the two wells and A is the tunneling split- 
ting originating from the overlap of the wave functions. 
The randomness of the glassy structure results in broad 
distributions for these two parameters. The phenomeno- 
logical tunneling model successfully describes most of the 
thermal, elastic and dielectric properties of glasses at low 
temperatures ||. However, experiments below 100 mK 
have revealed that deviations from the predictions of the 
tunneling model exist. Although a comprehensive theory 
is still missing, it seems that the interaction between the 
TSs gives rise to many of these deviations || . 

In recent experiments D. Rosenberg et al. investigated 
the AC dielectric constant during the first few sweeps of 
an electric bias field after the sample was cooled from 
about 1 K to about 10 mK |Q. The electric bias field 
was swept in triangular waves for two successive sweeps. 
The electric bias field amplitude of the second sweep was 
twice as large as the amplitude F\ of the first but the 
sweep rate dF/dt = Fi/tu with the period 2tn of the 
first triangular wave, was kept constant. At the begin- 
ning of the first sweep the dielectric constant suddenly 
decreased but then increased slowly with increasing elec- 
tric field. When the field reached the maximum and 
started to decrease again the dielectric constant suddenly 
increased and then decreased slowly with decreasing field 
to about its original value when the field reached F — 0. 
The solid line in Fig. [j] illustrates the dielectric constant 
during the bias field sweep. The slow increase/decrease 
of the dielectric constant while increasing/decreasing the 
bias field follows the dipole-gap behavior found in earlier 
experiments [Q. However, the sudden changes of the di- 
electric constant at F — and at F = F\ are remarkable 
new features. Typically the relative change of the dielec- 



w 
CO 



I 

8e*/e 

1 




ZFi 



electric field F 

FIG. 1: The change of the dielectric constant is plotted over 
the electric bias field. The arrows show the behavior in time. 
In a first sweep the field is swept in a triangle wave between 
zero and F\ . The full line illustrates the change of the dielec- 
tric constant during the first sweep. The dashed line represent 
the behavior during the second sweep where the maximum 
field is twice as large as in the first sweep but the sweep rate 
is the same. 



trie constant of these sudden changes was 5e'* /e' ~ 10 5 . 
It decreased inversely proportional to the temperature 
between T = 5 mK and T = 20 mK, where the effect 
could not be resolved anymore, be'* je' also increased 
with increasing sweep rate. The most remarkable fea- 
ture yet happened when the second sweep with twice 
the amplitude was performed. The dielectric constant 
smoothly increased with increasing field until F = F\. 
At this point the dielectric constant suddenly decreased 
to the value it had when this field was applied during the 
first sweep. Further increasing the field, the dielectric 
constant increased following the course which we would 
have extrapolated from the first sweep. The dashed line 
in Fig. [I] illustrates the dielectric constant during the 
second sweep. Reaching the maximum field 2F\ the di- 
electric constant suddenly increased and then decreased 
slowly with decreasing field. This behavior suggests that 
the sample remembers the field strength which was pre- 
viously applied to it. Only if the sample was left at zero 
bias field for several days or warmed above 1 K would 
the memory as described vanish. 
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An electric field F couples to the dipole moment p 
of a TS, thus, changing the asymmetry energy A(F) = 
A + pF. Since the distribution of asymmetries is flat 
over an energy range exceeding typical values of pF we 
expect in thermal equilibrium no changes of the dielec- 
tric constant with applied field. However, we expect the 
changing field to drive TSs out of equilibrium. The TSs, 
which dominantly contribute to the dielectric constant, 
have relaxation times r typically shorter than r = f s 
at the temperatures of interest H . Accordingly, isolated 
TSs driven out of equilibrium cannot be responsible for 
the observed long time effects. We propose that weak 
interactions between TSs cause the memory effect. 

Electrically allowed transitions lead to a dielectric re- 
sponse inversely proportional to the energy splitting and 
proportional to the occupation difference between the 
two states. Thus, small energy splittings result in large 
contributions to the dielectric constant. On the other 
hand, small energy splittings lead to almost identical 
occupation numbers for the two states so that their re- 
sponse is suppressed M. However, a field sweep drives 
the TSs out of equilibrium. If two TSs have the same 
energy splitting E\ — E2 for a specific field Fpp a weak 
interaction J12 <C £a,£«2 between these two TSs lifts 
their degeneracy resulting in a small splitting Ao p < J\2- 
We illustrate this 'avoided crossing' in Fig. || where we 
plot the excited energy levels of two TSs versus bias field. 
With increasing field the excitation energy of TS 1 (solid 
line) decreases until it reaches its minimum A01 and then 
increases again. Note that the energy splitting of a TS 
varies with the electric field since the field shifts its asym- 
metry. TS 2 (dashed line) behaves the same. The levels 
cross at F — Fpp where these two TSs form a resonant 
pair (RP). The TSs forming the relevant RPs have en- 
ergy splittings E ~ 2k^T and in thermal equilibrium 
the occupation difference of the split levels is approxi- 
mately proportional to Aop/(fcBT). Regarding the small 
number of RP (about one per 200 TSs), the excess dielec- 
tric response is negligible. However, we will show later 
that TSs with long relaxation times form RPs during the 
field sweep which are far from thermal equilibrium and 
provide an excess dielectric response of the order of the 
observed effect. 

In order to understand how RPs yield memory for the 
previously applied bias field consider two TSs forming a 
RP where TS 1 starts in the ground state and TS 2 in the 
excited state. An adiabatic field sweep from F < Fpp to 
F > Fpp results in a flip-flop leaving TS 1 excited and 
TS 2 in its ground state (compare Fig. |J) . For an ensem- 
ble of such TSs the subensemble of TSs 1 interchanges 
the occupation number of its excited state with that of 
the ensemble of TSs 2. We further assume that the relax- 
ation time of TS 1 is short, so that it adjusts its occupa- 
tion number (of its excited state) according to thermal 
equilibrium to its momentary energy splitting, and the 
relaxation time of TS 2 is long, so its occupation num- 
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FIG. 2: The spectrum of two tunneling systems is plotted 
versus bias field. At F — Frp these tunneling systems form 
a resonant pair. 

ber stays unaltered during the experiment. The flip-flop 
exchanges the occupation numbers and TS 2 will store 
the equilibrium occupation number to field F = Fpp. 
When the field reaches F = Frp the second time, after 
it was swept to its maximum value and then decreased, 
TS 1 has again an equilibrium occupation number to the 
applied field F = Fpp. Now both TSs are in thermal 
equilibrium. All further flip-flops store the equilibrium 
occupation number into TS 2 and the situation stays un- 
changed. Only during the very first formation of a RP 
one TS is in thermal equilibrium. This produces the ob- 
served memory effect. 

In the following we quantify the above picture and es- 
timate the temperature and sweep rate dependences. We 
start by calculating the resonant contribution || to the 
dielectric response function of a RP 

x(i-i') = ([P(t),P(t')])Q(t-t'). 

Here (• • • ) indicates the thermal and ensemble averages 
and [•, •] the commutator. O is the Heaviside step func- 
tion and P = p\o ZI+P20 ' Z 2 the total dipole moment of the 
two TSs. The absolute value of P depends on the magni- 
tude and the orientation of the individual dipoles pi with 
respect to the external measuring field. The Hamiltonian 
of a coupled pair of TSs is given by 

A 01 Ax A 02 A 2 

n = —^-(Jxl + —&zl H &x2 + —&z2 + J\20 z \Oz2 

with the Pauli matrices a a j describing the two state vari- 
ables of the two TSs. Considering only lowest order con- 
tributions in the interaction we obtain for the excess re- 
sponse xrp °f a RP 

Xb.p= (pi^+P2^) (^e) (n + -n_)2sin(^) 

(1) 

with the occupation differences n+ = tanh(E+ / (2k bT)) 
between the ground state and the second excited state 
of the RP and n_ = tanh(E_ / (2k B T)) between the 
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ground state and the first excited state. The energies 
E± = (E 1 +E 2 )/2±E P with E p -- 



A 2 p + {E X - E 2 )y4 
are the excitation energies of the RP and Ao p = 
(A i/-Bi)(Ao2/-B 2 )Ji2 (see Fig. |). Note that especially 
asymmetric TSs, Aj ~ Ei, forming RPs, have an excess 
dielectric response whereas isolated TS only contribute 
substantially to the resonant contribution of the dielec- 
tric constant when they are symmetric, A <C Ao. 

At first we estimate the response of RPs in thermal 
equilibrium. The standard tunneling model assumes a 
distribution function P(A, Ao) = Po/Ao leading to a flat 
distribution of energy splittings. One further assumes 
a maximum splitting E rnax ~ 1 K. The number #n of 
relevant RPs, formed by two thermal TSs, Ei ~ 2k B T, 
where the coupling fulfills 2Ao p > \Ei — Ej\, is j^n = 
{P a 2k B T) ■ (P 2A 0p )/(P Q E max ). 

A common estimate for the mean coupling between 
TSs are J ~ 10 ^mK jjj. Due to the l/r 3 -dependence 
of the interaction on the distance r between the TSs the 
mean interaction is proportional to the density of TSs 
and thus, the mean interaction of thermal TSs is J(T) « 
100 k B fiK at T = 10 mK with E max = 1 K. 

Assuming further Aj ~ Aoi for the TSs forming RPs 
and E p Ei, E 2 we roughly estimate the real part of the 
dielectric constant (given as the complex Fourier trans- 
form of the response function) for RPs (Aq p ~ E p ) at 
T = 10 mK to be 



l(eq) 
Xrp 



4 ■ 10" 4 ■ P p 2 

/(eg) 



with the mean dipole moment p. Xrp is negligible com- 
pared to the dielectric response, about ~ Pop 2 , of isolated 
TSs in the tunneling model. 

During a field sweep most TSs will not be in thermal 
equilibrium and we have to investigate the occupation 
numbers of two TSs which form a RP at a field Fpp dur- 
ing the sweep. For fields F ^ Fpp the coupling between 
the two TSs is irrelevant. Assuming a typical dipole mo- 
ment of 1 Debye, typical bias fields F\ — 133 kV/m 
change the asymmetry by an amount corresponding to 
a temperature $/k B — pF\jk B ~ 30 mK. With sweep 
times tp = 10 s the time, during which two TSs fullfil 
the RP condition 2Ao p > \Ei — Ej\ during the sweep, 
is 33 ms at T = 10 mK. This time is short compared 
to relaxation times of thermal TSs. Accordingly we ap- 
proximate the occupation difference n + — n_ of the split 
levels by n+ — ri_ ~ n\ — n 2 . Thereby, 



= e Jo 



ds ^ {s] {^(0) 



+ j ds^(s)eIo ds '^ s 'h a nh(E l (s)/(2k B T))} 



i(t) 



are the occupation number of isolated TS during an adi- 
abatic field sweep with constant sweep rate and with 
the initial occupation difference 71^(0). For an isolated 
TS typical experimental bias field sweeps are adiabatic 



if A > y l (pFi/k B )(h/t R ) ~ 400 nK with Planck's 
constant h. The flip-flop dynamics of a RP only oc- 
curs when the formation of the RP is adiabatic as well 
and, thus, the two TSs, which form a RP, must fullfil 

A 0p > J {pFi/k B )(h/tp) ~ 400 nK. The relaxation rate 
7i(s) of a TS is given by the one phonon rate 

7<(») = r-\s) = 70 A^( s )coth (J^) (2) 

with 70 = (1/cf +2/cj )B 2 /(2irph 4 ) where c/, t is the longi- 
tudinal and transversal speed of sound, p is the mass den- 
sity of the glass and B the strain coupling constant [[|. 
The time dependence is given by the time variation of the 
asymmetry energy due to the bias field. We neglected 
additional relaxation processes which might emerge from 
the interaction between the TSs R] since their relevance is 
still a matter of debate. Including them will change quan- 
titative results but not the qualitative picture evolved. 

Statistically a RP A with tunneling splittings Aoia 
and A02A has a partner RP B whose tunneling splittings 
Aois and A02B are chosen to match Aois = A02A and 
A02B = Aoia- The initial asymmetries have to obey 
AiB = Aia and A 2 b = A 2 a in order that both RP are 
formed at F = Fpp. According to Eq. (|l|) the responses 
of RPs A and B differ only by their occupation numbers 
and we obtain for the combined response of both RPs 

x' — Xo{(niA - n 2B ) + (niB - n 2A )} , 

where all factors except the occupation numbers are in- 
cluded in xo and the occupation numbers are considered 
at F = Fpp when the RPs are formed. Each of the occu- 
pation differences uia — n 2B and n± B — n 2A only depend 
on a single tunneling clement. Thus, we can discuss these 
two RPs as if we had two RPs each formed by two TSs 
with the same energy and the same tunneling splitting. 

To obtain a substantial dielectric response the RP must 
be formed with an occupation difference close to one. If 
both tunneling elements are the same the only difference 
between TS 1 and TS 2 (compare Fig. ||) results from 
the bias field sweep; the energy splitting of TS 1 is only 
decreasing until the RP is formed whereas the energy 
splitting of TS 2 goes through a minimum. Accordingly, 
we expect a maximum value for the occupation difference 
of TSs whose relaxation time is similar to the time the 
energy splitting of TS 2 is smaller then k-sT . Note that 
the relevant RPs are formed by thermal TSs E ~ 2k-oT. 
In Fig. [5] we plotted n\A — n 2B as a function of Ao for 
various energy splittings E. 

Numerically we find a maximum for 711,4 — n 2B for a 
tunneling splitting AJ with 

- . 1 - * x 2kpT 
t(A*(E),E) = -;f—-tR 

PF max 

The full line in Fig. ^ represents n\A — n 2B for an energy 
splitting E = 2k B T and the dotted and dashed line have 
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FIG. 3: The occupation difference of uia — «2B is plotted 
versus the tunneling splitting for a temperature T = 10 mK 
and pFup/ks = 3 K. 



E > 2k B T and E < 2k B T respectively. RPs with E 
2k B T contribute maximally and we define 



A* = A* (E = 2k B T) 



lpF„ 



1 



t R 2 7o T 2 



(3) 



Note that the adiabatic assumption for the RP restricts 
the tunneling splitting of the slow relaxing TS to values 
A > 90 yK where we used 70 = 10 8 K 3 s"\ E max = 1 K 
and J(T) = 100 k B fjK at T = 10 mK. Rosenberg et al. 
jjj find a minimal tunneling splitting Aomin = 250 fiK 
justifying the adiabatic approximation. 

To estimate the sudden increase of the dielectric con- 
stant 5s'* js' when the field reaches the maximum we have 
to count all RPs formed by one TS with Ao ~ AJ and a 
second TS with A < AJ. The values for tl\a — n 2B and 
n\B — n 2 A are taken from Fig. ^| and we obtain for the 
change of the dielectric constant 



5s'* 2 k B T ( A* 

—T = NrpXrp * 1-6 ■ Pop 2 In — -2 



(4) 



with the average dielectric response of a contributing RP 

2p 2 1, v 
XRP l(^T)2 [niA - n2B)ma * 

and the number of contributing RPs 



p3K B T p3A* /2 pE+J pA* 

N RP = dE dA I dE' I dA' 

Jk B T JA*/2 



Pn 



E-J JA 0mi „A-o (PoE max )A' 



PoE„ 



A Or, 



where we approximated P(E,Aq) ~ Po/Ao since Ao <IC 
E 1 for the relevant tunneling splittings. We estimated 
the number of TSs as PoE max . Eq. (^) predicts a lin- 
ear increase up to a maximum and then a linear decrease 
of 5s'* js' with increasing temperature. Eq. (^) further 
predicts a logarithmic dependence on the sweep rate. 



Both predictions are in partial agreement with experi- 
ments and we refer to Rosenberg et al. |l| for a detailed 
comparison. 

Qualitatively the dynamics of driven RPs explain the 
experiments except for the initial sudden decrease. The 
response due to the driven RPs is always positive. Thus, 
our model predicts a sudden increase in s' at the be- 
ginning of the first sweep and another such an increase 
when the field starts decreasing. The experimental ob- 
served initial decrease might be explained by considering 
the AC measuring field. Typically, the energy shift by the 
AC field pF AC is of the order of 0.24 - 2.4 mK k B |J by 
assuming dipole moments of 1 Debye. At temperatures 
between 5 — 20 mK one might accordingly expect the 
AC measuring field to drive the TSs substantially. Thus, 
the DC field provides only a 'second' sweep for the RPs 
formed near zero field. As soon as the DC field exceeds 
the AC field all RPs are formed the first time yielding a 
sudden decrease at the beginning of the DC field sweep 
as experimental seen. The typical energy splitting of RPs 
is small compared to the energy shift induced by the ex- 
perimental AC field. Therefore the AC field influences 
the flip-flop dynamics of the RPs. A future investigation 
should take into account the AC field beyond linear re- 
sponse but this is beyond the scope of the present letter. 

In conclusion, we have shown that a first sweep of 
an electric bias field prepares resonant pairs of tunnel- 
ing systems, which are temporarily formed during the 
sweep, in metastable states. Resonant pairs in these 
metastable states significantly contribute to the dielec- 
tric constant leading to a higher dielectric constant in 
subsequent sweeps, and thus yielding a memory for pre- 
viously applied fields. 
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There is an ensemble of TSs for each set of parameters 

Ao and A. By using the phrase 'occupation number' for 

a single TS we refer to the ensemble of systems with the 

same parameters. 

At the temperatures T < 20 mk relaxational contributions 
to the dielectric constant are expected to be negligible. 



